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Abstract 

We establish an integral representations of a right inverses of the Askey- Wilson finite dif- 
ference operator in an L? space weighted by the weight function of the continuous g-Jacobi 
polynomials. We characterize the eigenvalues of this integral operator and prove a g-analog of 
the expansion of e*^^ in Jacobi polynomials of argument x. We also outline a general procedure 
of finding integral representations for inverses of linear operators. 
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1 Introduction. 



The Askey- Wilson divided difference operator P^, [Q is defined by 

(1.1) iVJ)ix)= ■ a; = cos^, 
where 

(1.2) {5,g)ie^')=g{q'/'e^')-g{q-'/'e^'). 

Observe that z(g^/^ — g~^/^)sin^ which appears in the denominator of (1.1) is 5qX, x being the 



identity map: x ^-^ x evaluated at x. Magnus ^Tj showed how the Askey- Wilson operator arises 
naturally from divided difference operators. 

It is easy to that 

g"/2 _ q-n/2 
T^qTn{x) = ^1/2 _ g-1/2 

where T„(x) and Un{x) are Chebyshev polynomials of the first and second kinds, respectively. 
Therefore T>q maps a polynomial of degree n to a polynomial of degree n — 1. As such T>g resembles 
the differential operator. It was observed in ^ and |]l5l that one can construct integral operators 
which are a right inverse to Vg on certain weighted spaces. In |TH| Ismail and Zhang diagonalized 



the right inverse to Vq on weighted L2 space on [—1, 1] with Jacobi weights (1 — x)"(l + x)^ or 
g-ultraspherical weights 

1 1 - 2(2x2 - l)g" + g2n 

Recall the notations 

n 

(1.3) (a; g)o := 1, (a; g)„ := J]^ (1 - ag^), n = 1, 2, . . . , or 00, 

m 

(1.4) {ai,...,am; q)n= Yiia-k] q)n, 

k=l 

for g-shifted factorials. We shall normally drop "; g" from the shifted factorials in (1.3) and (1.4) 
when this does not lead to any confusion. Thus for the purpose of this paper we will use 

(1.5) {a)n := ia;q)n 



2 



A basic hypergeometric series is 



(1.6) r4s 



1 ai,. 












[ h,. 


■,bs 



q, z] = r(psiai, . . . ,ar;bi, . . . , b^; q, z) 



(^_l)"g"("-l)/2 



{q,bi,...,bs;q)n 
A g-analog of Jacobi polynomials is given by 



l+s—r 



(1.7) Pt^\x-q) 



{q, -q; q)r 



-4*^3 



q"+\ -q^+\ -q 



where x = cos6. The P„'s are called the continuous g- Jacobi polynomials. In what follows we 
assume 



;i.8) e'' = X 



+ y^^^^, e-'^ = x- V^^^, 



where the sign of the square root in (1.8) is taken so that — 1 x as x — > oo in the complex 
plane. The normalization in (1.7) was introduced by Rahman in The original normalization 
used by Askey and Wilson in ^ is P^'^\x\q)^ where 

(1.9) P^^)(x; q) = ^ ^q—Pi-^^\x\q'). 

(-?; q)n 

Both {PjC'^\x; q)} and {P^"'^)(x|g)} are called continuous g- Jacobi polynomials because they are 
orthogonal with respect to an absolutely continuous measure. The orthogonality relation is 

(1.10) Pi^'^\x\q)P!C'^\x\q)w^Ax\q) dx = h';^'^\q)Sm,n. 



The weight function Wa,i3{x\q) and the normalization constants h^^'^^q) are given by jT^, (7.5.28), 
(7.5.30), (7.5.31)]. 



;i.ll) Wa,f3{cos9\q'^ 
;i.l2) hl^'^\q) 



^ ^ga+l/2gj6l^ ga+l/2g-i6»^ _g/3+l/2gj6l^ _g/3+l/2g-ie . 
2vr(l - ga+/3+l)(g(-+/3+2)/2^ g(a+/3+3)/2. 



(?, g-^+l, _g{"+/3+l)/2^ _g(a+/3+2)/2. 



a+1 „/3+l _^(a+/3+3)/2. ^n(2a+l)/2 



(g"+\g'^+S-g 



(1 _ ^2n+«+/3+l)(^^^ qa+l3+l^ _g(a+/3+l)/2 . 
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Note that (7.5.31) in |10| contains a misprint where "g"(2«+i)/4" ^^ie right-hand side should read 

)5^n(2a+l)/2)! 

The Askey- Wilson operator acts on continuous g-Jacobi polynomials in a very natural way. It's 
action is given by, |jlO|, (7.7.7)], 



r)„-n+(2o+5)/4/'i „a+/3+n+l\ 

.ox r) p(^'P)(r\n) - ^ " {^-<1 ) p(a+i,m)/ I X 



Following [|r^ we define densely on L2\w a^fi[x\q)\ by 

oo 9^-n+(2a+5)/4/i ^Q!+/3+n+l\ 

U^J A. (^^^(„+/3+l)/2)(l + ^(a+/3+2)/2)(l ^^1^^' 

if 

oo 

(1.15) f-Y.fnP^n^'\A<l)- 

n=0 

Clearly Vq as defined by (1.14) and (1.15) maps a dense subset of L'^a^fi{x\q)\ into 
L2[wQ,+i,/3+i(a;|g)]. We are interested in finding a formal inverse to "Dg, that is we seek a linear 
operator Ta,/?;^ which maps L2[wa+i,/3+i(a;|g)] into L2[ifa,/3(x|g)] such that "Dg T^a,p;q is the identity 
map on the range of P^. It is clear from (1.14) and (1.15) that we may require Tq,^^.^ to satisfy 

°° n _oU_o(a+/3+l)/2.ol/2N n-(2«+l)/4 

(1.16) (r„,,, ,)(x) ~ E ^ 2ri-a-^'/^'4 (-1^) 

n=0 ^l-^ y 

for 

oo 

(1.17) g{x)^Y.9nP't^^^^^^\A^)- 

n=0 

It is easy to find a representation of Ta^p-q as an integral 

operator. We use the orthogonality relation (1.10) to write Qn as 

(?(x)i/;„+i,^+i(x|g)p("+i'^+i)(x|g) dx/h^^-''^^~''\q) 

and formally interchange summation and integration in (1.16). The result is the formal definition 

(1.18) (T„,^.g^)(x) = J^^Ko,,/3-g{x,y) g{y) Wa+i,/3+i{y\q)dy, 
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and the kernel Ka,f3;q{x,y) is given by 

The purpose of this work is to study the spectral properties of the integral operator (1.18). It 
is worth noting that T^^^.g is linear but is not normal and a spectral theory of such operators is not 
readily available. Our main result is Theorem 3.1 which characterizes the eigenvalues of Ta^p-q as 
zeros of a certain transcendental function. In order to prove Theorem 3.1 we proved several auxil- 
iary results which may be of interest by themselves. First in §2 we solve the connection coefficient 
problem of expressing Pj{^'^\x) in terms of {Pj°'^^''^^^\x)} . The solution of this connection coeffi- 
cient problem is then used to expand Wa+i,i3+i{x\q)P^'^^^'^'^^\x\q) in terms of {Wa,f3{x\q)Pl:"'^\x)}. 
In Section 3 we used the latter expansion to find a tridiagonal matrix representation of Ta^p-q. In 
Section 3 we also find the eigenvalues and eigenfunctions g{x\X) so that 

Ta,l3;q 9 = ■ 

The eigenvalues are multiples of the reciprocals of the zeros of a transcendental function X^^{^\l/x) 
defined in (5.13). Such a function is a g-analog of the confiuent hypergeometric function iFi. The 
eigenfunctions are shown to have the orthogonal expansion 

oo 

(1.20) Y.^^{\\q)PlC^^\x-q) 

1 

such that 

oo 

(1.21) 5:e'^)(g)|a„(A|g)r < oo. 

1 

We may normalize the eigenfunctions by choosing ai(A|g) = 1. With this normalization we prove , 
in Section 3, that a„(A|g) is a polynomial of degree n. The a^'s are q-analogs of polynomials studied 
by Walter Gautschi [|lT] and Jet Wimp ^3]. An explicit formula for a„(A|g) is given in Section 4, 



see (3.11) and (4.1). The large n asymptotics of a„(A|g) in different parts of the complex plane are 
found in Section 5 and they are used to characterize the A's for which (1.21) holds. 

In Section 5 we note that when a and (3 are complex conjugates and are not real then the 
polynomials i~^an{ix\q) are real orthogonal polynomials. They are orthogonal with respect to a 
discrete measure supported at the zeros of a transcendental function which we denoted by Fiirj, p; q). 
The function Fiijq^p^q) is a g-analog of the regular Coulomb wave function Fiirj^p), |l|, Chapter 
14]. Our analysis implies that the functions FlItj, p; q) have only real and simple zeros. Neither the 
functions F^ij]^ p; q) nor their zeros seem to have been studied before this work. 
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In Section 6 we prove that the eigenfunctions are constant multiphes of the g-exponential function 
Sq{x] —i,b), where 

(1.22) S,ix; a, h) := ^ ^{aq^'-^^^'e^' , aq^'-^y^e^^Xb". 

n=0 W/" 



The function £g was introduced in |15|. Since the eigenfunctions also have the orthogonal expansion 
(1.20) we obtain an identity valid on the spectrum of Ta^p-q. This identity is (6.13) and is a g-analog 
of 



(1.23) 







n + p + 1 
2n + a + (3 + 2 



2ty] P(°'^)(x) 



-1 < X < 1, 



see (10.20.4) in In Section 6 we use properties of basic hypergeometric functions to show that 
the above mentioned g-identity holds also off the spectrum of Ta^p-q. In Section 7 we give a second 



proof of (6.13) using a technique similar to what was used in to prove the same result for the 
continuous g-ultraspherical polynomials. We also include in §7 a formal approach to finding the 
spectrum of certain integral operators of the type considered in this paper. In §8 we include some 
remarks on asymptotic results of Schwartz Dickinson, Pollack and Wannier 0] and general 
remarks on this work. 



In many of our calculations we found it advantageous to follow ||T0[. The only disadvantage is 
that we have to introduce some additional relations. Recall that the Askey- Wilson polynomials are 
defined by [0, (7.5.2)] 



'1.24) j9„(x; a, 6, c, (i|g) = (a6, ac, Ofi; g)„a "4^3 



q ahcdq 



ab, ac, ad 



Their orthogonality relation is ^ (7.5.15), (7.5.16)], 
/■I /^(x;l,-l,gl/^-gl/2) 



;i.25) 



h{x; a, b, c, d) 



-Pn{x; a, 6, c, d)pm{x; a, 6, c, d) 



dx 



x^ 



q;q 



^ ^ ^1 ^{1 — abcd/q){q,ab,ac,ad,bc,bd,cd)n ^ 
' ' ' {1 — abcdq^^~^){abcd/ q)n 



with 



^1.26) /i(cos^^; ai, 02, 03, 04) := (0^6*^,0^6 



4 

id „ „-ie\ 

00 ; 
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and 

(1.27) K,{a, b, c, d\q) = 2n {abed) oo/[{q, db, ac, ad, be, bd, cd)oo]- 
We shall also use the notation 

(1.28) /i(cos 6; a) := (ae'^, ae'"'')^. 
Note that 

(1.29) h{cose;l,-l,q'/^q-'/^) - {e^'',e-^''U. 

Note also that the continuous g-Jacobi polynomials correspond to the identification of parameters 

(1.30) a = ?(2«+l)/^ 6 = 5(2a+3)/4^ c=-?(W/4^ ^ _^{2/3+3)/4_ 

In fact the polynomials Pjl"'^\x\q) of (1.7) have the alternate representation, 

^^a+l.^^^ / q-n qU+a+P+l ^(2a+l)/4gje ^(2a+l)/4g-ie 



;i.31) Pt''\x\q)^ 



q;q 
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2 Connection Coefficients 



In this section we derive a g-analogue of the formula 

(2.1) (1 - a:^)Pin^''^'Hx) = , 4(n + a)(n + /j) ( 

4n(a-/3) T^ia^B)^^ 4n(n + 1 



(2ri + a; + /?)(2n + a + /3 + 2) " (2n + a + /? + l)(2n + a + /? + 2) 

Formula (2.1) is (1) in §37 of |19[. The Jacobi polynomials {P^"'^^(x)} are orthogonal with respect 
to (1 - x)"(l + x)^ and (2.1) is essentially the expansion of (1 - x)°+i(l + x)^+ipi"t^'^^^Ha;) 
in terms of {(1 — x)°(l + x)^Pj"'^\x)}'^Q. The g-analog of this question is to expand 
Wa+i,(3+i{x] q)Pi"X^'^^^\x; q) in terms of {Wa,i3{x] q)Pj"''^\x; q)}j^Q. It is worth mentioning that 
the latter problem is equivalent to the expression of P^"'^-*(x|g) in terms of {Pj"^^''^~^^\x\q)}^^Q. 
This follows from the following known observation. If {pn{x', A)} are orthonormal with respect to a 
weight function w{x; A) then the connection coefficient formula 

n 

Pn{x; X) = Y1 Cn,j(A, yU)Pj(x; /i) 
j=0 

holds if and only if its dual, namely 

oo 

w{x] fJ.)pn{x; /i) ~ 51 Ci,n(A, fj.)w{x; X)pj{x; A), 

j=n 

holds. This latter fact follows from computing the Fourier coefficients of both sides. 
The main result of this section is the following theorem 

Theorem 2.1 We have 

(2.2) (1 - 2a;g"+i/2 + g2"+i)(l + 2xqf'+'/^ + q"'+')PtV'^^'\x; q) 

_ (1 + q^+P+^)(l + ^a+/3+n+l)(l + g"+'^)(l + g/5+n)(l _ g«+«)(l _ g/3+«) 

" n-'i 9') 



;i + g°+/3+"+l)(l + ga+/3+2n+l)(^ ^ gU^^^ _ gU^^^ _ ^a-P 
(^\ _ q2n+a+l3^(^Y — g2n+Q:+/3+2') 

:i + g")^(l + g"+^)^(l-g")(l-g"+^) 

(1 _ g2n+a+/3+l)(-]^ _ g2n+a+/3+2) ^ '^+1 \-^->^>- 



■qPplC^P\x-q) 
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For our purposes it is convenient to express (2.2) in the Askey-Wilson normalization P^"''^)(a;|g). 
The result is 

(2.3) (1 - 2xg"+i/2 ^ ^2a+i)(^ ^ 2xg^+i/2 ^ ^2/3+i)p^°+i./3+i)(a;|^2) 

- (g2n+a+/5; g)2 ^ ^"-1 ^^'^ ^ 

, (-g"+^+^ g)2(l + g"+/^+^"+^)(l - g^")(l - q-^) .^^n-ipj.,,).... 



(g2"+"+/3+l; g)2 



-r 



[x\q 



The rest of this section will be devoted to proving (2.2). Our proof is very technical and the 
reader who is more conceptually oriented is advised to turn to Section 3. 

Our proof uses the Sears transformation [0, (III. 15)] 



(2.4) 403 



q "',a,b,c 
d,ej 



q "',a,d/b,d/c 



{e,f)n \ d,aq'--/e,aq'-^/f 



where def = abcq^ Note that the parameters q '^,a and d remain invariant under the transfor- 
mation (2.4). 

We now proceed with the proof. We seek a connection coefficient formula of the type 

gi-", 6cg"+^ gVSg-ie 



(1 - 26cos0 + r)(l + 2ccos^ + c^)40, 

~h ''"^'y hqy\-cqy\-q 

where h := g""*"^/^, c := 

The orthogonality relation (1.23) gives 



—c(fl'^, —q 



AfeK(g^/^6, -c, -g^ 



/2 ^ (1 - hc){q, -he, -b^, c^)k 
(1 - bcq'^''){bc, -g, ~Cy/q, by/q)k 



/g2ie -2.e^ 



h{cose;^,bq,-cq,-y^) 



■ 493 



q ^,bcq'^,q^/'^e^^,q^^'^e 
6g^/^, —cq^/"^, —q 
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gi-", 6cg"+^ gi/^e*^ g^/^g-is 



g, g (i6l. 



Apply the Sears transformation (2.4) with invariant parameters g 6cg*^ and — g to the first 4^3 in 
the above equation. The result is 



(2.5) «:(gV2,6,_c,-gV2) 



(-g)^(l -&c)(g, -&c)fc 
(1 - 6cg2fc)(6c, -g)fc 



At 



' (g-^ bcq'')r q' (g^"", 6cg"+^), g^ 

^ „^ ^^1/2 _A^l/2^ 
r= 



{q, -Q, cq^^^, -bq^/^)r J^q {q, -q, bq^^^, -cq^^^). 



/■71 

Jo 



de 



h{cos9; g^+^/2, bq, —cq, — g^'+^Z^) 

The integral on the right-hand side of the formula (2.5) is 

27r(&cg'-+^+3)^ 

(g, 6g^+3/2, -cg^+3/2, -g'^+^+i, -6cg2, -6g-+3/2, cg-^+^/s)^' 

which can be written as 

.1/2 , ^ 1/2^ (l + M(l + M(l-g&^)(l-gC^) 

K[q , 0, — c, — g J — — — — 

{bcq)2{bcq^)r+s 



■{-bq'/',cq'/%.{bq'/\-cq'%i-q) 



r+S' 



Therefore (2.5) leads to 



(2.6) 



(l-feg)(l-fecg^) il-bc)iq,-bc)k 
[-60)2(1 - qb^){l - qc^) (1 - bcq^^){bc, -g)^ 



^ {q~\bcq\-bq^'\cq^l%. , 
.t^o (g,&cg3,cgV2,-6gi/2)^ ^3</>2 



gi-", 6cg"+\ -g^+i 



6cg 



r+3 



The above 3^2 can be summed by the g-analog of the Pfaff-Saalschiitz theorem, (11.12) in [|T^]. It's 
sum is 

(-6cg^g^+'-")„-l 
(-gi-", 6cg'-+3)„_i 
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which clearly vanishes if r < n — 2. Thus = Q ii k < n — 2. When k > n — 1, replace r by 
r + n — 1 in (2.6) and simplify the result to see that the right-hand side of (2.6) is 



(-?!-", hcq^ cq^l\ -bqy^)n-i {bcq-+^)n-i 



When k = n — 1 the above is 1. If A; = n, n + 1, the aforementioned 4^3 has only 1, 2 terms; 
respectively. After some simplification we find 

_ (1 - &gV^)(l + cgV^)(l + bcq'-){l + &cg")(l + g-)(l - b/c) 

_ (1 - 6gV2)(i + cgV2)(l - 6g-+V2)(i ^ cg"+V^)(l + g")(l + g"+^) 
^^•^^ (l-6cg2n)(l-6cg2n+i) ^^^^ ^' 

_ (1 - bq^'^){l + cg^/^)(l + 6g"-^/^)(l - cg"-^/^)(l + &cg"-^)(l + 6cg") 
^ ■ ^ ~ (l-6cg2-i)(l-6cg2-) ■ 

Using the dual relationships mentioned at the beginning of this section we see that Theorem 2.1 
is equivalent to the following theorem. 

Theorem 2.2 The connection coefficient formula 

„-n/2('i „a+/3+n+l\/i „Q!+/3+?i+2\ 

f2 10) P^^'^VxIo) = - - — - - P^^+^-'^+^^fxIo) 

y^.LV,) yxiq) ^_^(„+^+i)/2. ^i/2-)^^2 - g"+("+/3+l)/2) (1 _ gn+(a+/3+2)/2-) K-^IU) 



^(.+;3+2-n)/2(i _ ga+/^+n+l)(i ^ ^n+(a+„a+l)/2) _ g(«-/3)/2) ^^^^^^^^ 
+ (^_^(a+/3+l)/2. ^1/2)2(1 _gn+(a+/3)/2)(]^_^n+(a+/3+2)/2) "^"-1 V^l^^ 



^(3a+/J+4-n)/2(i _ ^a+n)(i _ 



(_^(a+/3+l)/2. 51/2)2(1 - 5»+(«+/3)/2)(l _ 5n+(a+/3+l)/2) 

holds. 



P^V'^^^ixlq), 
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3 An Eigenvalue Problem 



In this section we characterize the eigenvalue and the eigenfunction of the eigenvalue problem 

(3.1) {T^,p,,g){x)^\g{x). 

This characterization is stated as Theorem 3.1 at the end of the present section. 

It is tacitly assumed in (3.1) that g belongs to the domain of T^j^.g and \g belongs to its range. 
Now assume 

oo 

(3.2) g{x) g{x- \\q) ~ ^ an{\\q)Plt^^\x\q). 

n=l 

Since g e L2[wa,i3{x;q)] then (1.10) implies 

oo 

(3.3) E h(^'''\q)MX\q)\' < oo. 



n=0 



The condition (3.3) and the eigenvalue equation (3.1) will characterize the eigenvalues A and the 
eigenf unctions g. 

It is clear that 

oo 

(3.4) AEa„(A|5)p("'/^)(x|5) 

=1 

oo 1 

= Ea„(A|g) / Wa+i,p+i{x-q) K^"'^\x,t)Pt'^\t\q) dt 

n=l -^-^ 

^ (g,g-^^g^^^-g^,-g^))oo(l-g) ^-,,..,n 



n=l 



■2^an{\\q)2^ q {x\q) 

n=i fe=o ^ -g 2 )jt 

• j]w^+i,p+,{x-q)Pl'''''^'^'\t\q)Pt'\t\q) dt. 

The next step is to evaluate the integral on the extreme right-hand side of (3.4). This will lead to 
a three term recurrence relation satisfied by the a„'s. By (2.3), the integral on the right, denoted 
Ik,n with k replaced by A; — 1, is 

(3.5) 7fe,„ = W(t; q-/^+y\ g-/2+3/4^ _//2+l/4^ _//2+3/4) 
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a + P+3 . 



(1 - - q 2 

(1 - g 2 - g 2 



Now the orthogonality relation (1.10) is equivalent to 

w{t- q'-l^+^l\ q-/^+^l\ _g/3/2+i/4^ -qP'^+^'^)P^-^^\t\q)P^'''^\t\q)dt 

01+I3+2 a+13+3 

2TT{q 2 ,g 2 



a+/3+3 



(1 _ ^a+/3+l+2,-)(^«+/3+l^ _^(a+/3+l)/2. 

Substituting this in (3.5) we find that 

(3-6) /m = 7 IT-^T^ i+HI a+/3+2 , (1 + g 2 )(1 + ? 2 ) 

(g,g"+i,gP+i,-g 2 , -g 2 ;g), 



loo 



(l-g"+^+i) (g"+\g^+S-g^;g)n .p^+i) 

U 5 ^ J (?"+'^+\?, 2 ;g)„ 

(l_ga+fc)(l_g/3+fc) ^ 
(1 - g 2 +'=)(! - g 2 +'=) 

a + /3+l I , , , , a-t3 , 

(1 + g 2 — — g 2 ) /3-a+fc-l50z49z 

^ (l-g^+^)(l-g^+^) ^ ' 
(l-g^)(l-g^+^)g^^^ 



From (3.4) and (3.6) we have 

00 

(3.7) Ai:«n(A|g)P("''^)(a;|g) 



n=l 



(l-g)(l-g 2 )(l-g 2 ) .^vi + 

2(l-g"+i)(l-g.«+i)(l-g"+/3+3)^ ^ ^ ^ > 
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oo -, Q,+/3+i+2fc / a+/3+3 H±|±3 \ , ^ , , „^ 

fc=i J- 9 ^ 2 



(1 _ qa+(3+l+2k^ 



After some simplification we find 

oo 

(3.8) AEan(A|g)Pi"'^)(x|g) = Y.Pt^\A<D 



oo 



n=l fc=l 



2(1 - g 2 +'^)(1 - g 2 



_ (1 - - - ,i3+t+i),(3»+g+t+i)/2 

2(1 - g"+m+fc)(i _ g^^+fc)(i - g^^+fc) 

By equating the coefficients of P^'^\x\q) on botfi sides of (3.8) we establisfi the following three-term 
recurrence relation for a„'s 

MM f+i (1 - g)(l - g"+^+^)(l - g^+^+^) 

2(1 - g°+/3+l+fc)(l _ +^)(1 - 5^^+*=) 

(l-g)(l-g^)(l + g^+^) ^ 



2(l-g— +'=)(l-g 



g 2 afe(A|g) 



(l-g)(l-g"+^+^) ^ 



2(i-g^+fc)(l-g^^+fc) 
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The limiting case g — > 1 of the recursion relation (3.9) is 

2{a + l + k){(3 + l + k) 



(3.10) -Aafc(A) 



{a + P + 1 + k){a + p + 2 + 2k)2 
2{(3-a) 



afc+i(A) 



-afc(A) 



2{a + P + k) 



-afc_i(A). 



{a + p + 2k){a + P + 2 + 2k) {a + P + 2k - l){a + P + 2ky 

for A; > 0, which is (4.16) of [|l^ 

It is clear from (3.3) that ao(A|g) = and that ai(A|g) is arbitrary. It is also clear from (3.3) 
that ak{X\q) /ai{X\q) is a polynomial in A of degree k — 1. It is more convenient to renormalize 
ak{X\q) / ai{X\q) in terms of monic polynomials. Thus we set 



(3.11) ak+iiX\q) 



, ^+13+5 , 



.l)'^'5j.(^:^^)g-('=V4+("+/3/2+l)fc)_ 



{q»+\qf'+';q)k ^ ^'^l-q 

There is no loss of generality in taking 6o(2A/(l — q)) = 1. In terms of the 6„'s, (3.8) becomes 



(3.12) bk+iifx) = bkifi) 



1 - g 2 )(1 + g 2 



a/2+3/4+fc/2 



0(1 







_ g" + l + ^)(l — g/J+l + fc^g's+^ + l 



(1-g- 



'')(1 - g^^+*')2(l - 



where 



(3.13) /^=^^, &-i(/i)=0, 6o(/i) = l. 

In Section 5 we shall determine the large n behavior of the polynomials &„(x) and a„(A|g). These 
asymptotic results will be used to prove the following theorem. 

Theorem 3.1 The eigenvalue problem (3.1)-(3.3) has a countable infinite number of eigenvalues. 
The eigenvalues are (1 — g)/2 times the reciprocals of the roots of the transcendental equation 



(3.14) (-p"+3/2(i_^)^/2;p)^20i 



p 



;i - g)xpi/V2 



1/2 



-j9°+3/2(l - g)x/2 

Furthermore A = zs not an eigenvalue and the eigenspaces are one dimensional. 



Theorem 3.2 Any eigenfunction g{x; X\q) corresponding to an eigenvalue X is a constant multiple 
of£q{x; -i,X). 
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4 A q- Analog of Wimp's Polynomials. 



In this section we find an expficit solution of (3.12). 



Theorem 4.1 The polynomial generated by (3.12) and (3.13) are given by 

-i)v/V""' 



(4.1) 64/x) = y: 

j=0 



^p-l3-n-l^_p-a-n-l.p-^. 



403 



,J2n+a+l3+3-j „/3+l _„a+l 
f -I F -1 1/ -I 1/ 

a+13+2 n+l3+2-j a+n+2-j 



P^' ' ,p 



P" 



where 



(4.2) p := q 



1/2 



Proof. From (4.1) it is clear that bQ^fi) = 1 and that 6i(//) satisfies (3.12) when k = 0. Since the 
solution of the initial value problem (3.12)-(3.13) must be imique, all we need to do is to verify that 
the right-hand side of (4.1) satisfies (3.12). The actual process of this verification is rather long and 
tedious. 

First, we shall rewrite (4.1) in the form 

(4.3) bM = E E B^-k (-i)^-^V/V"-^ 

j=0 k=0 



where 
(4.4) Ak 



(p^+\ -P"^';p)k ^(n) _ (p-^-"-S -p-"-"-^p)fc 



{p,p- 



-2n-a-f3-2. 



;p)k 



Since Aq = Bq"'^ = 1, we then have 



n+l j 



(4.5) br.+^{^^) = + E E (-1)^+V/V"+^- 

j=l k=0 



n J 



bnin) - E E i-iy^'p^/'n-' + E E i-iy^'^p'/'n-^'-^ 



n+l j 



(n+l) 



j=l k=0 



j=0 k=0 
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where the last hne is obtained by separating /x" from the rest of the series on the right-hand side of 
(4.1). Our first aim is to bring the same factor of bn{(J^) as is shown in (3.12), so we make a further 
separation of the series on (4.5) and find that 



(4.6) bn+i{fi) = bM 
where 



n—lj+1 n j+1 

j=0 k=0 j=0 k=0 

[i- P A^+P )P W A, f-1 V+'^n^//"-^ 

The rest of the exercise is to show that c„(/x) is actually a multiple of 6„_i(/x), the same multiple 
as in (3.12). The coefficient of jjP- in Cn(/x) is 

which vanishes by the use of (4.4), verifying that Cn(/i) is a polynomial of degree n — 1 in We 
thus have 

r4 81 C fu) - ^ '-- - V ^i-lV+'^AuB^^ o0+2)/2„n-j-l , J / N 

where 

n-2 j+2 

(4.9) ci„(/.) := E E(-l)^'^'^^^fc<2V^'+^' V-^-^ 

j=0 fc=0 



n-1 i+2 

- E E(-i)''"''"'^^^iSilp^^'^'^/V"-^-^ 

j=Q k=Q 



Since 



j(n) 75(n+l) _ IP . P>j+l-k 



. ■^-2n-a-^-3^-|^ -p'"*^+^)(l _ p-2n-a-/3-4-j _ ^-n-^-2^-|^ -p'®"")(l _ p-4n-2a-2/3-5+j-fe^ 
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we find that the coefficient of /x" ^ ^ in is, for 

< j < n - 2, 

(4.11) (Lj^ — i:/'=«&"(-i)'"V'=-(i-/-')p-''-"-= 

. V 4,n _ „-4n-2a-2/3-5+j-fcw_-,\j+fe 0+2)/2_lP ,P)j+l-k 

Now we combine the second series in (4.11) with the coefficients of /j,"^'^'^ in the ffist series on 
the right-hand side of (4.8) which, by virtue of the identity 

(^l _|_ pa+/3+3+2n^^Q+n+l ^-^-n-2^2 _ p-4n-2a-2/3-5+j-fe'^ 



(1 _ pa+/3+2+2n) (1 _ pa+/3+4+2n) (^-2n-a-/3-2. p).^^_^ (p-2n-a-/3-4. p).^^_^ 

(^l _ p?+l-fcjp-/3-"-2 
= -(1 -p«+/3+2+2n)(p-2„-a-/3-3.p)^.^3_^' 

results in the series 

(1 _p-/^-n-l)(l J o(n-l) , . ,/2 

(l_p2„+a+/3+2)(p-2n-a-/3-3.p)3 A. ^ j-fc ^ ) V ■ 

Adding this to the first series in (4.11) we find, after a straightforward calculation, that the coefficient 
of < j < n - 2 in c„(/i) of (4.8) is 

(-\ „2a+2n+2\/i „2,(3+2n+2\„a+/3+2n+2 j 
y^-"-^) (^1 _p2n+a+/3+3j^l _p2n+a+/3+lj(-l _p2n+a+/3+2^2 Z_^^'=-^j-fc ^ > ^ ' 

Finally, collecting the j = n — 1 term from the series in (4.8) and (4.9) we find that the constant 
term in (4.8) is given by 



f4 13l P A^+P )P A g(n) /_iY«+fc-l (n+l)/2 



n+l 



k=0 



say. By (4.4), we get 

(p;p),+l(p-2"-"-/3-4;p)„+2 



(4-14) /n = .1^.. . ' J-.-.-. (-l)n-V+l)/2 
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'n~a-l3-3\ 



„-n-l „n+Q:+/3+4 „/3+l 



a+2 



P,P 



(^1 _ pQ+/3+2+2n^ 



a-2/3-3n-6 



„— n „n+a+/3+3 „/3+l 

P 1 P -if 1 P 



a+2 



P,P 



To simplify the expression on the right side of (4.14) we first denote 



F ^ 1/ ) 1 F 

„a+/3+2 /3+2 _„«+2 
F -I F -I f 



p,p 



and then use the contiguous relation of Askey and Wilson, see |T0|, Ex 7.5] 



B C 

(4.15) 0„+i = -^0n - ^0n-l, 

where 



(4.16) A = p°+/5-3"+4(^]^ _p"+a+/3+3^(^]^ _ p2n+a+/3+2^^-|^ _ ^„+a+/3+2^ 



(l-p"+^+2)^l^pn+a+2^)^ 



n+a+2\ 



5 = -C - A + p"+/3-3n+4(p2„+a+/3+2)^(^ _ p/3+l)(l + 



Substituting (4.15) in (4.14)we find after some simplification that the coefficients of 0„ cancel out, 
so that 



(4.17) fn = -ji-ir-V""^'^^' 



(p-/3-n-2^_p-»-n-2.p)^^^ 
(p,p-2n-a-/3-4.p)^_^^ 



■>n-l 



'1 _ p2n+a+l3+3'j(^l _ p2n+a+/3+l^ (-j^ _ p2n+a+f3+2y 



k=0 



Combining (4.12) and (4.17) we find that c„(/i) is the same as the second term on the right side of 
(3.12) (with k replaced by n). This completes the proof of Theorem 4.1. 
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5 Properties of {bn{x)}. 

In this section we derive asymptotic formulas for the polynomials {bn{x)} in different parts of the 
complex x-plane. We also investigate a closely related set of orthogonal polynomials and record 
their associated continued J-fraction. 

Recall that we normalized {bn{x)} of (4.1) by 
(5.1) boix) := 1. 

To exhibit the dependence of bn{x) on the parameters a and P we shall use the notation 6^"'^)(x) 
instead of bn{x). Our first result concerns the limiting behavior of {b^^'^\x)}. 



Theorem 5.1 The limiting relation 

(p'^+i,-p°+3/2x;p), 



(5.2) Wm x''b^^^^\l/x) 



holds uniformly on compact subsets of the complex x-plane. 



Proof. We may use Tannery's theorem (the discrete version of the Lebesgue bounded conver- 
gence theorem) to let n —>■ oo in (4.1) after multiplying it by Thus the sequence {x~"'6"'^(x)} 
will have a finite limit if the series 



^ (-a;)-V'/2 

7—^ 3</^2 



0, 



p,p 



converges. Therefore 

(5 3) hm (X) - T y (^-^^^"^-^°""^)^ 

(5.3) „limx 6„ (x)-l^^ ^^.^^^ p-H-^^-p), ^' 

if the right-hand side exists. In the above sum interchange the j and k sums and replace j hjj + k 
to see that the right-hand side of (5.3) is 



-p' 



'-/2 



The j sum is (p^^^/x; p)oo by Euler's sum ||10|, (II. 2)]. This shows that 



(5.4) lim bl^'f^ (x) = (p'/^x; p)oo20i 



^pa+l pf3+l 



P 



a+13+2 



p.p^^"^ /x 
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uniformly on compact subsets of the open disc {a; : |a;| < p^^^}. On the other hand from Theorem 
8.1 we know that 

hm x"6l"'^)(l/x) 

exists uniformly on compact subsets of the complex plane and is an entire function of x. The Heine 
transformation, |IU| , (III.l)] 



(5.5) 201 (a, b] c; q, z) = ^f' ^^^°° 20i(c/&, z; az; q, b), 



(c, z,q)c 



implies 



_pa+l pf3+l 



p 



a+13+2 



P,P 



,/3+l 



Therefore (5.2) holds in the interior of {x : |a;| = p^/^} and analytic continuation establishes the 
validity of (5.2) on compact subsets of the complex plane. This completes the proof. 

We next determine the asymptotic behavior of h^^'^^[x) at x = and in : < |x| < p^/^}. 



Theorem 5.2 We have for a ^ P, 
(5.6) 6^^)(0) ^Cp^'/'w" 



as n — >• cxo, 



where C is a nonzero constant and \u\ < 1. 



Proof. Clearly 



(5.7) 6i"'^)(0) 



(p-^~"~\-p'"'"'^p)n 
(p,p-2«-°-/5-2;p)„ 



P 1 P 1 P 1 P 

F ? F 1 P 



-i)V^^ 



Ismail and Wilson [14| proved that if < 1 then 



(5.8) 403 



q ^^ahcdq^ ^,az,a/z 
ah, ac, ad 



a \ " {az, bz, cz, dz; q)c 
z) {z^, ab, ac, ad; q)c 



as n ^ oo. Now apply (5.8) with a = ip^H^+P)/"^^ b = -ipi+(°+/3)/2, c = -ip^-i^-'^)/'^ , 

d = ,pl-{"-/3)/2^ 

^p(a-/3)/2 i{a> (3 
_ipiP-a)/2 if ^ > a. 
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Therefore 



{az, bz, cz, dz; q)c 



Zj {Z\ p°+/3+2 , pP+2 ^ -pa+2 . ^ ^ 



which imphes (5.6). 

Corollary 5.3 The values of C and u in (5.6) are given by 



(5.9) u 



(l+p<.-3)(pa + /3 + 2.p)^ tj a ^ P 



c 



(i+p/3-«)(p"+f+2;p)^ 13 > a- 
The only case left now is to determine the large n behavior of 6^"'^-'(0) on the zeros of 



Q+1 „l/2 



X 



,/3+l 



This is not straightforward and requires some preliminary results. 



Theorem 5.4 The function 



(5.11) Yt''\x) = {-xy 



a+2+k 13+2+k 



■2(Pl 



P 



a+l3+2k+A 



P 



1/2 



X 



satisfies the three term recurrence relation (3.12), with p = q^^'^ . 

To prove Theorem 5.4 we used MACSYMA to first find the multiple of the 20i then verified 
that the Y^"''^^ of (5.11) indeed satisfies (3.12) by equating coefficients of powers of 1/x. Theorem 



5.4 is also a limiting case of a result of Gupta, Ismail and Masson [|T^, as will be explained in §8. 
It readily follows from Theorem 5.4 that 



Xl-^^\x) = i-x)-'{p'/'/x;p)^2^^ 



_pa+2+v pl3+2+u 



P 



a+l3+2v+A 



P 



1/2 



X 



satisfies the three term recurrence relation 



(5.12) 



(1 _ gV+(a+/3+3)/2-)(^]_ _ gU+{a+l3+4) /2y(^l _ ^,.+(a+/3+5)/2) 



^ I l-^ g Jl-^ + g )_ {u+a+3/2)/2 

(1 _ qU+ia+l3+2)/2^(^l _ qU+{a+l3+4)/2^ 
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The Heine transformation (5.5) yields the alternate representation 

(5.13) = (-.y '^f 



_pa+u+5/2 



According to Theorem 4.5 of [15 



(5.14) aa+i ■ ■ ■ a+n-iXSf (x) = h^t^^f'^^\x)X^:^^\x) + er^''^^'^^'H^)4-f ^(a:) 

where is the coefficient of Xyj^\{x) in (5.12), see (4.28) in [|T5[ . 

Theorem 5.5 Let ^ be a zero of F(x) of (5.10). The large n behavior ofh^^'^\^) is 



(5.15) e'^)(o 



(^(a+/3+3)/2^g(a+/3+5)/2.g)^ (g(a+/3+4)/2. g)2^ X^^'^^ (0 ' 



Proof. It is clear from (5.10) and (5.13) that F{x) = if and only if X^{^\x) = 0. The recurrence 
relation (5.14) shows that if Xjf'^\x) and xl^f\x) vanish at x = C then xl'^n\0 = for all n, 
n = 2, 3 But it is obvious that 



[X) ^ {—X) , as n —> oo. 



Thus Xj^"'^)(x) and X^"'f''(x) have no common zeros. Now (5.14) implies 

" (g("+/3+3)/2,g("+/3+5)/2;g)^(g("+/3+4)/2.^)2^X("'^)(0' 

and we have estabhshed (5.15). 

We are now in a position to prove Theorem 3.1. 

Proof of Theorem 3.1. From (1.12) and (3.11) it follows that 

(5.17) hi-f,\q) K^,{\\q)\' = 0(g-"('^+^"+^^-^/^)/>(-^)(2Ag^/V(l - q))]'). 

Now (5.2), (5.6) , (5.9) and (5.15) show that ET ht^K^l) MMq)? converges if and only if 
2g-'^/^A/(l — g) is a zero ofxi"/^(x). The eigenspaces are one dimensional since the eigenfunc- 
tion with an eigenvalue A must be given by 

oo 

(5.18) g{x\\)=Y.aM<l)Pt^\x\q)- 
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Next we prove that the polynomials {i "'b^°'''^\ix)} are orthogonal on a bounded countable set 
when a and P are not real and are complex conjugates. 

Set 

(5.19) 4"'^)(a;) :=r"6(f'^)(«a;). 

The three term recurrence relation (3.12) leads to the following three term recurrence relation for 
the s„'s. 



(5.20) stfi\x) 



X + 



^g{a-l3)/A _ ^(/3-a)/4^^(2n+a+/3)/4 



^(1 _ gn+l+{a+(3)/2-^^l _ ^„+2+{q+/3)/2) 

_ (1 - g"+"+l)(l - gn+/3+l-)gn+l+(a+/3)/2 ^^^^ 

~ (1 _ gn+l+(a+/3)/2)^^n+2+(a+/3)/2.^1/2)g'^"-l ^^f' 

We also have the initial conditions 

s[''^^\x) :=0. 



X) 



i-'^\x) 



(5.21) s. 
When 

(5.22) a = p, Im a 7^ 0, Re a > -1, 

then the coefficient of sl^'^\x) in (5.19) is real for n > and the coefficient of sl^li\x) is negative 
for n > 0. Thus the s„'s are orthogonal with respect to a positive measure, say dip. The coefficients 
in the recurrence relation (5.19) are bounded. Thus we can apply Markov's theorem |21[], namely 

#(t) 



(5.23) lim(.i"'^)(x))74"'^Hx) 



Im X 7^ 0, 



' — oo X — t 

where (s^"'^^(x))* is a solution to (5.20) satisfying the initial conditions 

(5.24) {s^^'^\x)r := 0, (4"'''^(x))* := 1. 
It is easy to see that 

(5.25) (sl-'f'Kx))* = st\'''^'\x). 

Therefore (5.2), (5.18) and (5.23) give 

diP{t) (^"+^^+2.^)2 



(5.26) 



X 



t (l-p/3+l)(l_^pa+3/2/3.) 



201 (p"+^ - V^Vx; ip"+^/Vx; p, 

201 (p"+^, —ip^/'^x; ip'^'^^/'^ / x; p, pf^~^^) 
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Recall that the Coulomb wave function FL{ri,p) are defined in terms of a confiuent hypergeo- 
metric function as [0] 

(5.27) F,(ry, p) := ^'e^^^'^-^Yi^^^^ + 1 - ^r^; 2L + 2; 2^p) 



As g — > 1 the polynomials (1 — g)"s^"'^''(x/ (1 — q)) tend to the Wimp polynomials, |2^. The right 
hand side of (5.2) in the case of Wimp's polynomials is Fl{t],p). This suggests defining a g-analog 
of Fl{ti,p) by 



(5.28) Fi(r/,p;g):=(V/V;g)oo20i 



_„L+jj7+l ,.L-iri+l 
H 1 H 



q,iq^''^p 



where L and t] are real parameters. Observe that the iterate of the Heine transformation [|T^, (III. 3)] 
(5.29) 2<l)i{a,b]c;q,z) = (^^^/^' (-^^ 5. ^^^^ _ 



shows that -F/;,(?7, p; g) is real when p is real, as in the case of Fl^t], p). 
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6 An Expansion Formula. 



The purpose of this section is to give a direct proof of the eigenfunction expansion (6.13). 
Set 

oo 

(6.1) £g{x; a,r) = ^ amPmix; b, -c, -c^). 

m=0 

In order to compute am we need to evaluate the integrals 

(6.2) Jm(a; r) := J ^ w{x; b, b^/q, -c, -Cy/q)pm{x; b, b^, -c, -c^/q)8q{x] a, r)dx, 

when a = —i. We shall keep the parameter a in (6.2) free till the 

end then we specialize the result by choosing a = —i. It is clear that 

(6.3) Jm{a; r) = ^ — — — Jm,„(a, b, c), 

n=0 W' 1)n 

where 

w(x; 6, 6^, -c, -c^)pmix; b, fe^g, -c, -cyg)— ^ — j+^dx. 
-1 /i(x; ag 2 ] 

Formulas (6.3.2) and (6.3.9) in |jlO| imply 

w{x;a,/3,-f,6)— — — dx 
-1 h{x]f) 

^ 27r(a^, (3g, 7^, /ff, a(3-f6f/g; g)oo 
(g, a/3, 07, a(5, a/, /?7, /?/, 7^, 7/, 5/, g)^ 

■8W^7(fi'V?; 5f/a, 5(/7, ^f/ /; g, a(3-f6f/g). 
Using the 4(^)3 representation of p„ in (6.4) we obtain 

27r(a6g(i-")/2^ afegi-'^/^ -acg^i-'^)/^ -acgi-'^/^. 



(6.5) Im,n{a,b,c) 



(g, fe^^g, — &c, —bcy/q, —bcy/q, —bcq, (^^-i abq^^^'^^l'^; q)c 
(ga^62c2g"+^;g)oo 



E 



{abq^+^'l'^, -acq^'^+^'il'^, -acq^+^'l'^, a^q^-^] q). 



j=o {q,b^c^q''+\abq(^-''y^)j 
■8iy7(a'g""; ag-^-("-i)/V^ ag^^^V^, -ag(^-")/Vc, -ag-"/Vc, g""; g, b^c\^+''+^). 
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We now apply Watson's transformation formula which expresses a terminating very well-poised s4>7 
as a multiple of a terminating balanced 403, |]10|, (III. 17)]. Thus 



We then apply the Sears transformation (2.4) with invariant parameters g — ag "^^/c, g ""^^^"^/c^. 
After some simplification we obtain 

(6.6) 8W^7(a'g-"; ag-^-("-^)/V6, ag~"/V6, -ag(^-")/Vc, -ag-"/Vc, g""; g, fe'c^g^+^+i) 
(a^gi-", c2gi/2, -g-"/6c, -g-"+i/V6c)„ 



(-acg(i-")/2, -acgi-"/2, -a6g(i-")/2, -a6gi-"/2). 



(_5cg-+V2^a5g(i-n)/2). 



(a6g("+i)/2,-6cgV2), 



•4<P3 



g-", g-"-V&^c2, -ag-"/Vc, -g-"/Vc 
-g-"/6c, -g-"--''+i/V6c, g-^+^/Vc^ 



9,? 



The substitution of the right-hand side of (6.6) for the gW^y in (6.5) and some simplification lead to 



Im,n{a,b,c) = K{b,c)- 



^g^/^, —bcq^^"^, —hcq)n 
{qb'^c^)n 



{-a/cTq 



{q,b^c^q''+\-bcqy^), ^ 



g-", g-"-V&^c2, -aq-'^/^/c, -g-"/V 
-q-^'/bc, -g-"--''+i/V6c, g-"+i/Vc2 



ac 



where 

(6.7) /t(6,c) 



2Ti{qb\^) 



(g, 62gi/2^ —bcq^/"^^ —bcq^/'^, —bcq, c^q^/'^)^ 
27r(6cgV2,6cg)oo 



(g, 62gi/2^ —bcq^^'^, c^q^^'^)oo 
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Replace the 403 by its series definition with summation index k. Then interchange the j and k sums 
to obtain 



" {q-\ -aq-^/yc, -g-^'^/ac, g'^/b^^ ^ 



fc=0 



3<P2 



,2„2„n+l-fc 



The 302 can now be summed by the g-analog of the Pfaff-Saalschiitz theorem |T^, (11.12)]. It's sum 
is 



which vanishes for all A;, < < n if m > n. Thus we get 



..ON J i (g-, -g-+VV&c)^ , (c^gV^-6cgV^-&cg)n , , ..2/2 



(-6cgi/2,g-™-™/52c2)^ 



g'"-", g-'^-"/62c2, -ag-"/Vc, -g~''/'^/ac 
-g-"/6c, -g-"+i/V6c, g-^+^/Vc^ 



The relationship (6.3) and the observation lm,n{,0'^ 6, c) = if n < m show that Jm(a; r) is given by 



^m?/A m 00 /„„m/2\n 

E /'t^^^ n+^(a, 6, C). 



Applying (6.8) we find after some simplification 

(6.9) .Ua,r) = .(^ c),"V.(_,f,)™(£VV^^^zM 

(go C )2m 



~ (c2g'"+l/2^ _5cg"^+l/2^ -bcg"'+^] 

h (g,62c2g2-+i)„ 



ar 



-m/2 



g-'^, g-2™-"/62c2, _ag-{«+m)/2/c, _g-{«+™)/2/ac 
-g-^-'^/fec, -g-"-™+i/V6c, ^-"-"^+1/2/^2 
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The 403 on the right hand side of (6.9) has a quadratic transformation. By (3.10.13) of [IC], the 
aforementioned 403 is equal to 



_q-n-m 1^^^ _^(-n-m+l/2)/2 j ^(-n-m+l/2)/2 



(gl/2, -g-"-™/6c, g(-n-™+l/2)/2/c, _g(-n-m+l/2)/2/c; ^ 



g-"/2^ _5cg™+("+l)/2, cg(''"+l/2)/2^ _cg(™+V2)/2 



^^^m+1/2^ -cg('"+i)/Va, -acg 



(m+l)/2 



by reversing the sum in the first 4^3. After some straightforward manipulations we establish 

(6.10) Jm{a,r} = K{b,c) . q ' 

[bcqy^, bcq; g)„ 

og^ (-cg(-+i)/Va,-acg(-+^)/^;gV^)„ / ar _^/,_y,Y 
h {q'l\-q^l^:q^l^)n \ 



g-"/2^ _5cg™+{«+l)/2^ ^^(m+l/2)/2^ _Cg(™+l/2)/2 j 



Finally we apply Sears transformation (2.4) to the 4^3 in (6.10) with invariant parameters q 
^g,(m+i/2)/2 hccf^^^l'^ . This enables us to cast (6.10) in the form 



(6.11) J^(a,r) = K{h,c) ^ q ' 

[hcq^l^hcq; q)m 

g (-agV^-gV7a;gV^). ^^^^„ 



n=0 



(gV2,_gl/2.^1/2)^ 



• 4<P3 



g-"/2^ ^^(m+l/2)/2^ _5g(m+l/2)/2 

5cg-+V2^ _ag(-"+l/2)/2^ _^(-n+l/2)/2/^ 



^1/2^^1/2 



It is evident from (6.11) that Jm{a\ r) is a double series. When = —1 we have been able to reduce 
the right-hand side of (6.11) to a single series. To see this, replace the 4^3 in (6.11) by 

its defining series then interchange the sums. The result is 



{bcq^/^,bcq;q)m {-ir;q) 
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bcq 



m+l/2 



We have tried to express Jm(a; r) of (6.11) as a single sum for general a but this does not seem to 
be possible except when a = ±i. 

Now the orthogonality relation (1.10) gives for the case a = —i 



(6.12) K{b, c)am = Jm{-i; r) 



;i_i,2^2^2m)^^2^2^^2^1/2^_^^.^) 



m i-2m 



(1 — b'^c'^)(q, c^q^/"^, —bcq; q). 
Thus we established the expansion formula 



(6.13) 8q{x] -i, r) = amPm{x; b, bq^^'^, -c, -cg^^^), 

m=0 

with the a^'s given by 



,1/2 
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7 A Formal Approach 



We now formalize the procedure followed in Section 3 and used earlier in ||T^. Let S G and assume 
that for every A & S, the sequence of polynomials are orthogonal with respect to a 

measure with a nontrivial absolutely continuous component. By A+1 we mean (1 + ai, . . . , 1 + a^) 
ii A = (ai, . . . , Ofe). We will assume that A + 1 E S whenever A E S. Let the orthogonality relation 
of the PnS be 



(7.1) 



Pn{x; A)pm{x; A)diJ,{x; A) = h^iA) 5m,n- 



Assume that T> is an operator defined on polynomials by linearity and by its action on the basis 
{pn{x; A)} via 

(7.2) Vp^{x;A)=UA)pn-i{x;A + l). 

Furthermore assume that the support of fi'{x; A) = tell is the same for aW A E S and that we 
know the connection coefficients in 

n 

(7.3) pnix; A) = Y^ Cnj Pjix; A + 1). 

j=0 



dx 



Therefore 
(7.4) c„ 



1 



^"'^ hj{A + I) J ^oo 

The formula dual to (7.3) is 



Pn{x] A) pj{x; A + 1) dfi{x; A + 1). 



(7.5) pn{x; A + 1) i^'ix; A + l)= £ wp„(x; A) //'(x; A), 

holding on the interior of the support of yu'(x; A). 

We now wish to describe the spectrum of a formal inverse to T). Note that we can define V 
densely on L'^{dfi{x; A)) by (7.2) provided that the polynomials Pn{x; A) are dense in L'^{d^{x; A)) 
for all A E S. This suggests that we define via 

oo oo 

(7.6) J2 cinPn{x- A + 1):=Y. anPn+i{x- A)/in+i{A). 



n=Q 



n=0 



This motivates the definition 



(7.7) (T^/)(x) := 



fit) 



E 



Pn+l{x] A) Pn{t] A + 1) 



^0 in+i{A)K{A+l) 



d^{t]A + l) 
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if / e L^{dij{x;A + l)). 

The next step is to consider the eigenvalue problem 

oo 

(7.8) TAg = Xg, g{x)^Y. ^niX; A) A). 

n=0 

In order for (7.8) to hold it is necessary that g lies in the range of T4, hence (7.6) shows that 
ao{X;A) = 0. Therefore (7.7) and (7.8) yield 

f a.(A; A) p„(x; A) = ± ^^±ii^ £ ^f^'f H p.(t; A)p„(t; A + l)d^^{t; A + 1). 

n=l n=l fc=„/lnl^+ij -^-oo 

Thus we have established 

(7.9) A Y: «n(A; A) p„(x; A) = ^ T\\\ E a,{X, A). 

n=l n=l Sn+ll^J 

Now (7.9) implies the recurrence relation 

00 

(7.10) XUA)an{X;A)= ^ 

{X,A). 

k=n—l 

Observe that (7.10) transformed the eigenvalue problem (7.8) to the discrete eigenvalue problem 

(7.10) . When ^ = for k < n — r for a fixed r then (7.10) is the eigenvalue equation of a matrix 
with at most r + 1 nonzero entries in each row. The cases analyzed in and in this paper are 
the cases when 

(7.11) Cn^k = for k < n — 2, n = 2,3, ■ ■ ■ . 

Note that c„,„ 7^ 0. When (7.11) holds then (7.10) reduces to 

(7.12) XUA)a^{X;A) (A; A) + (A; A) + 

Cn+l,n— l^n+l 

{X-A). 

For example in the case when the p„'s are the ultraspherical polynomials C^{x) we have |T^, , §144], 



(7.13) ^C:{x) = 2uC:ll{x), 2(n + z/)C:(x) = 2z/[Cr^(x)-C:^2^(x)]. 
Therefore 

^n{j^) = 2Z/, Cn,n = -Cn,n-2 = 1^ / {v + u), Cn,n-1 = 0, 
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and (7.12) becomes 

a„_i(A;z/) a„+i(A;z/) 



(7.14) 2Xan{X]u) 



{v + n-l) [v + n + iy 
which is (2.11) in 

The procedure just outhned is very formal but can be justified if both p„(x; A) and p„(a;; A + 1) 
are dense in L'^{dfi{x; A)) fl L'^{dfi{x] A + 1)). 

Another approach to the same problem is to think of T4 as a right inverse to V. In other words 
VTa is the restriction of the identity operator to the range of T>. Thus (7.8) is equivalent to 

(7.15) g{x;X) = XVg{x;X). 

Now the use of the orthogonal expansion of g and formulas (7.2) and (7.15) implies that 
XC,n+i{A) a„+i(A, A) is the projection of Y^^Li afc(A; A) pk{x] A) on the space spanned by Pn{x; A+1). 
Therefore (7.3) implies (7.12). 

In [|r^ it was observed that the eigenfunction expansion Y.kLi '^fc(A; A) Pk{x; A) can be extended 
to values of A off the discrete spectrum of the operator under consideration. This is also the case 
with the expansion formula (6.13). We now attempt to find such an expansion in general. 

We seek functions {Fn{X; A)} such that the function£'yi(x; A), 

00 

(7.16) Ea{x; A) := Ffc(A; A) p,(a;; A), 

k=0 

satisfies 

(7.17) XVEa{x;X) = Ea{x;X). 

Observe that (7.16) reminds us of (7.15), hence under the assumption (7.11), (7.3) and (7.17) imply 
that the F^s must satisfy a recursion relation similar to (7.12), that is 

(7.18) XUA)Fn{X;A) = 

Cn— l,n— 1 Fn—\ Cn,n— 1 Fn {X-A). 

In order to maintain a parallel course with the results in [|I^ and with the notation of Section 5 
we will renormalize the a„'s and -F„'s in order to put (7.12) in monic form and change (7.18) to a 
recursion with the coefficient of y-n-i equal to unity. Keeping in mind that ao(a;; A) = and that 
ai{x\ A) is a multiplicative constant which we may take to be unity, we set 

(7.19) a„(A;A)= n^^^^w"~^&n-i(AM;A)n> 1, ai(x; A) = 6o(x; A) = 1. 

j=0 9+2,i 
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and 



(7.20) Fn{X; A) = H t^^j^ Gn-i{Xu- A),n>0 Fo(A; A) = G.i{Xu; A). 

j=o 

Here m is a free normalization factor at our disposal and may depend on A. Thus 

(7.21) A6„(A; A) = 6„+i(A; A) + 5„(A)6„(A; A) + C„(A)6„_i(A; A), 
and 

(7.22) AG„(A; A) = C„+i(A)G„+i(A; A) + i?„(/l)G„(A; A) + G„_i(A; A), 
hold with 

(7.23) i?n(A):=p±^, a(A) = "';"?;"7;; . 

We are seeking a solution to (7.18) that makes (7.16) converge on sets of A's containing the 
spectrum. Since the p„'s are given we need to choose the F„'s to be as small as possible, that is 
choose Fn to be the minimal solution, if it exists. Recall that a solution Wn of (7.18) is minimal 
if Wn = o{vn) where f„ is any other linearly independent solution of the same recurrence relation. 



IGj . It is clear that a minimal solution of (7.18) exists if and only if (7.22) has a minimal solution. 



The minimal solution may change form in different regions of the parameter or variable space, [p!2 
When Bn{A) —>■ and Cn{A) —>■ as n oo then (7.22) has a minimal solution, see Theorem 4.55 
in 



le 



In many cases we encounter a fortutious situation where we can choose u in (7.23) such that 

(7.24) Bn{A) = Bo{A + n), Cn{A) = Co{n + A). 
Therefore 

(7.25) G„(A; A) = Go(A; A + n) = G_i(A; A + n+1). 



When (7.24) holds Theorem 4.5 of comes in handy. This latter theorem is 
Theorem 7.1 Let /(x; u + A) he a multi-parameter family of functions satisfying 
(7.26) G,+^/(x; p + A + l) = {A,+ax + S,+^)/(x; u + A) ± f{x; u + A-l), 
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and let {fn,u+A{x)} be a sequence of polynomials defined by 



(7.27) fo^u+A{x) = 1, fi^u+A{x) = Ay^AX + By^A, 



(7.28) ,!/+a(3;) — (^n+i/+A3; + -Bn+!y+A)/ra,j^+A(3;) i Cn+i^+A-l/n.-l,j^+A(2;)- 



Then 



(7.29) C^+aC;.+a+i ■ ■ ■ C';.+A+n-i/(a;; z/+A+n) = fn,u+A{x)f{x, u+A)±fn-i,u+A+i{x)f{x; u+A-1). 
Theorem 7.1 establishes 

(7.30) C_i(A) C-i{A + 1) ■ • • C_i(A + n - l)Go(A; A + n) 



where we used G_i(A; A) = Go{\; A — 1), see (7.25). 

Now assume in addition to (7.26) that Bn{A) —>■ and Cn{A) 0, hence the minimal solu- 
tion to (7.22) exists. Let {Gn{)\] A)} be the minimal solution to (7.22). According to Pincherle's 
theorem, |jT6|, the continued J-fraction associated with (7.21) converges to a constant multiple of 
Go{X] A)/G^i{\; A). Therefore the eigenvalues of the infinite tridiagonal matrix associated with 
(7.21) are the zeros of G-i{\;A) which are not zeros of Go{\;A). If G-i{\;A) = then (7.30) 
indicates that the series in (7.16) is a multiple of the series (7.8) and the multiplier does not depend 
on X but may depend on A. This explains the relationship between the expansion representing the 
eigenfunction g in (7.8) and the expansion in (7.16) which is expected to be valid for a range of A 
wider than the spectrum of T>. 

Finally we apply the preceding outline to the case of continuous g-Jacobi polynomials and give 
another proof of (6.13)-(6.14). In the case under consideration 



bn{X; A)Go{X; A) + 6„_i(A; A + l)G'_i(A; A) 



(7.31) en(«,/3): 




(l-g)(-g("+/3+i)/2;gi/2)2 



(7.32) c, 



(_g(a+/3+l)/2. gl/2)2(l - gn+(a+/3+l)/2) (^^ - gn+{"+/3+2)/2) ' 



(7.33) c, 




■n,n—l ■ 



(_g(a+/3+l)/2. gl/2)2(l - gn+(a+/3)/2 ) ]^ _ ^n+(a+/3+2)/2) 
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and 



(7.34) Cn^n-2 '■- 



(_g(a+/3+l)/2. gV2)2(l - g'»+(«+/3)/2) (1 _ gn+(a+/3+l)/2-) " 

With the choice 

(7.35) = Y 

we find 

(7.36) 5.(a, /3) = Boia + n, P + n) ^ - '^^Ul^' LJ^J '"''^-''"^'' 

(1 - q 2 +'*)(! - q 2 +") 

and 

(1 _ g"+l+")(l - ^/3+l+nN^n+^4^+l 

(7.37) C„(a,/3) = Co(a + n,/3 + n) = .+^+1 .+g+2, «+^+3 , 

(1 - g 2 +")(!_ 2 +^y{l-q 2 +") 

Therefore (5.13) yields 

(7.38) ^^(A;^,/?) -G'o(A;a; + n,/? + n) = (-A)-("+^)/2^^"'^)(A). 

When A = {a,P) and Pn{x; A) arc the continuous g-Jacobi polynomials P^"'^^(a;|g) we will denote 
Ea{x; A) by Ea,i3{x; A). Now with V = Vg formula (7.16) becomes 

(7.39) E„A^; A) ^ g ^^jf;;.';;- ^^i;^^^ (r^)""G»-> (^;«,;3)p<°.'')(.i,). 

We now find another solution to (7.17) and prove a uniqueness theorem for solutions of (7.17). We 
then equate Ea,p{x]\) and the second solutions to (7.17) and establish (6.13). 

Lemma 7.2 Assume that f{x) is an entire function of the complex variable x. If 

(7.40) V,f{x)^'^f{xl 

then f{x) is unique up to a multiplicative function ofy. 



36 



Lemma 7.2 is essentially Lemma 3.5 in |15| 
A calculation gives 



(7.41) VySg{x]a,b) = — £g{x;a,b). 



Therefore Lemma 7.2 implies 

Theorem 7.3 /// satisfies the assumptions in Lemma 7.2 then 

(7.42) f{x)=w{y)S,{x--t,y/2). 

Theorem 7.4 The function [— 2Ag^''^/(l — q)Y°'~^^^^'^ Ea^/six; A) does not depend on a or (3. 
Proof. In general we have 

oo 

(7.43) Ea{x; A) = ^ F„(A; A) p„(x; A) 

n=0 

oo 

A + l) + 

n=0 
oo 

= ^Pnix; A + l)[c„,„F„(A; A) + c„+i,„F„+i(A; A) + c„+2,ni^n+2(A; A)] 

n=0 
oo 

= J2pn{x;A + l)\^n+iiA)F„.+i{X;A) 



n=0 

oo n 



\u Pn{x; A + 1) u--' n z^T^ UM) G„(Aw; A) 

n=0 j=0^j+l[^) 



= A ^ Cn,nPnix; A + 1) n ^^Gn-i{Xu; A + 1). 

n=0 j=OSj+U^j 

In the case of continuous g-Jacobi polynomials the last equation gives 



(7.44) E,,^(x;A) = A^ 



(_g{«+/3+l)/2. gl/2)2(l - gn+(a+/3+l)/2)(l _ ^n+(a+/3+2)/2) 

("^ 



? ^.fa+/3+l)/2.gl/2)2„ l^l-gj W-1 ^ _ ^ , « + i, P + i j ^„ l^igj. 
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Therefore 

The above functional equation can be put in the form 

[-2Agi/V(l - A) = [-2Agi/V(l - g)]("+^+^)/2E«+2,/5+2(x; A), 

and we have 

(7.45) [-2Agi/V(l - A) 

= hm [-2Agi/V(l-g)]'^'"-^("+^)/'i?2„.+a,2™+/3(a:;A). 

Now substitute the right-hand sides of (7.38) and (7.39) for Gn and E^^p in the right-hand side of 

(7.45) to get 

oo 

(7.46) [-2Ag^/V(l - g)]'"^+("+'^)/'i?2„^+.,2r.+/5(x; A) ^ A ^ g"'/^ (-A)-"g-""/^Pi"'^)(x|g). 

n=0 

But (1.24) and (1.31) imply 

„na/2 „na/2 

p(-'P){x\q) ^ ^-^p„(x;0,0,0,0|g) = -^if„(x|g), 

where {if„(x|g)}^ are the continuous g-Hermite polynomials. Thus the limit on the right-hand side 
of (7.45) exists and we have established 

(7.47) [-2Ag^/V(l - g)]("+^)/^K,;3(a:; A) = ^ f . H^Q)- 

n=0 W) 9Jn 

This proves Theorem 7.4. 
Corollary 7.5 We have 

(7.48) [-2Ag^/V(l - q)]^''^^^'^E^^p{x- A) = (A'^; g^)oo^,(x; -i, i/X). 



Corollary 7.5 follows from |15| where Ismail and Zhang proved that the right-hand sides of (7.47) 
and (7.48) are equal. 
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Theorem 7.6 The expansion ofSq{x; —i, r) in a continuous q-Jacobi series is given by ( 6.13) where 
6 = g(2a+i)/4 a„(/c = g(2/3+i)/4. 



Proof. Prom Theorem 7.3 and Corollary 7.4 we see that the right-hand side of (6.13) is 
w{r)£q{x] —i,r). Furthermore w{r) does not depend on a or /3 since neither Sg{x; —i,r) nor the 
right-hand side of (6.13) depend on a or (3. Now w can be found by letting a and (3 tend to oo then 
use (7.48). 
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8 Remarks. 



In 1940 Schwartz published an interesting paper ||20| containing the following result 



Theorem 8.1 Let {Pn^fi^)} be a family of monic polynomials generated by 
(8.1) po,uix) = 1, pi^u{x) = X + B^, 



(8.2) (x) + Cn+u Pn-l,u{x). 



If both 



(8.3) ^\Bn+u-CL\ < OO and ^\Cn+u\ < 



oo 



n=0 



n=0 



hold, then x" p„^;^(a + l/x) converges on compact subsets of the complex plane to an entire function. 



It is clear from (8.1) and (8.2) that {pn^u{x))* = Pn~i,u+i{x), hence the continued J-fraction as- 
sociated with (8.1) and (8.2) converges to a meromorphic function of l/x and the convergence is 
uniform on compact subsets of the complex plane which neither contain the origin nor contain poles 
of the limiting function. Schwartz illustrated his theory by applying it to the Lommel polynomials 
and he mentioned their orthogonality relation. 



Somehow Schwartz's interesting paper was not noticed and neither his results were quoted 
nor his paper was cited in the standard modern references on orthogonal polynomials 
and continued fractions, p2|, |1^. Many of Schwartz's results were later rediscovered by others. 
Dickinson, Pollack and Wannier rediscovered the special case Bn+u = of Schwartz's theorem. 



It is worth noting that if Bn+u = then a theorem of Van Vleck, Theorem 4.55 in []TB|, states that 
Cn+iy suffices to establish the uniform convergence of the continued J-fraction associated with 
(8.1) and (8.2) to a meromorphic function. The convergence being uniform on compact subsets of 
the complex plane which neither contain the origin nor contain poles of the limiting function. 



In |]T2[ it was proved that 



(8.4) 



X) 



xABc) (Ag", Cg", Dg^/A, Dg"/^, Dq'^/C), 



■3<P2 



Aq"", Bq"", Cq"" D 

^^2n^ -g"/x, xABC 
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satisfies the three term recurrence relation 



(8.5) Zn+l{x) = (X - an) Zn{x) - bn Zn-l{x), 

with 

(Ra^ n ._ D (1 - Dg'^mi - DqyB){l - Dq^C) 

^ ^ ABC ^ (l-L>g2n-l)(l_L)52n-2) 

D (l-Ag'^-i)(l-fig"-i)(l-Cg"-i) 



+ 



ABC (1 - L'g2"-i) (1 _ £)g2n) 



and 



(8.7) 6„ - ^?"-')(l - i^?'^-')(l - Cq^-'] 



(1 - L)g"-7A)(1 - Dq"-^/B){l - Dq''-^/C) 

■ (l-Dg2«-l)(l-Dg2n-2)2^1_2)^2n-3) " 

We next identify (3.12) as a hmiting case A — > oo of (8.5) with B — —C. When A — > oo, it is easy 
to see that 



(8.8) an 



(1 - Dg2n-2) (1 _ ^^2n) ' 
L»g2n-3 _ 52^2n-2) _ ^2^2n-2/52^ 



52 (1 - Dg2„-i) ^1 _ nq^n-2)2 (1 _ nq^n-3) ■ 

We replace 5 by g^/^ then identify the parameters A, B, C, D as 
(8.9) B = = -C, = g2+(a+/3)/2^ A ^ 00. 

It is not difficult to see that if Zn{x) satisfies 

Zn+l{x) = {X + a'^) Zn{x) - h'^ Zn-l{x), 

with 

g(«-l)/2(X + qn+(a+f3+3)/2-^ _ ^(/3-a)/2^ 



(8.10) a; 



(1 _ qn+l+(a+f3)/2'^ (^l _ ^n+2+(o+/3)/2) ' 
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^" ^ (1 _ gn+(a+/3+l)/2) (^l _ ^n+l+(a+/3)/2-)2 (^l _ ^n+(a+/3+3)/2^ " 

It then follows that 
(8.11) := g(2-+5W4 

satisfies (3.12) with /i replaced by x. This relationship between solutions of (8.8) and (3.12) enables 
us to take advantage of the detailed study of solutions of (8.5) contained in |12[. For example one 
can obtain the minimal solution to (3.12) by inserting the values of A, B, C, D of (8.8) into X^^^ of 
, which remains a minimal solution. This gives an alternate derivation of the form of Yj^°''^^ (x) 
of (5.11). 
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